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Motivated by recently discovered unusual properties of bulk nematic elastomers, we study a 
phase diagram of liquid-crystalline polymerized phantom membranes, focusing on in-plane nematic 
order. We predict that such membranes should generically exhibit five phases, distinguished by their 
conformational and in-plane orientational properties, namely isotropic-crumpled, nematic-crumpled, 
isotropic- flat, nematic-fiat and nematic-tubule phases. In the nematic-tubule phase, the membrane 
is extended along the direction of spontaneous nematic order and is crumpled in the other. The 
associated spontaneous symmetries breaking guarantees that the nematic-tubule is characterized by 
a conformational-orientational soft (Goldstone) mode and the concomitant vanishing of the in-plane 
shear modulus. We show that long-range orientational order of the nematic-tubule is maintained 
even in the presence of harmonic thermal fluctuations. However, it is likely that tubule's elastic 
properties are qualitatively modified by these fluctuations, that can be studied using a nonlinear 
elastic theory for the nematic tubule phase that we derive at the end of this paper. 
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I. INTRODUCTION 

Fluctuating poiymerized[| membranesQ have at- 
tracted considerable interests in the past two decades. 
Probably the most striking property that distinguishes 
them from their one-dimensional polymer analogue and 
most other two-dimensional systems, is that polymerized 
membranes admit a low-temperature flat phase 0, 0, 
characterized by a long-range order in membrane's nor- 
mal. This phase is separated from a high-temperature 
crumpled phase, with randomly-oriented normals, by a 
thermodynamically-sharp transition [(J 0,0- At low tem- 
peratures, strong interplay between thermal fluctuations 
and elastic nonlinearities infinitely enhances the effective 
bending rigidity, which in turn stabilizes the flat phase 
against these very fluctuations y, 0, la, u|- The resulting 
flat phase exhibits rather unusual length-scale dependent 
elasticity, nonlinear response to external stress and a uni- 
versal negative Poisson ratio. 0,0 

More recently, stimulated by a number of possi- 
ble anisotropic realizations, the role of elastic mem- 
brane anisotropy was considered by Radzihovsky and 
Toner 10]. It was discovered that while permanent in- 
plane anisotropy is unimportant at long scales within the 
crumpled and flat phases, it qualitatively modifies the 
global phase diagram of anisotropic polymerized mem- 
branes, leading to an intermediate (in its properties and 
location in the phase diagram) "tubule" phase. In the 
tubule phase the membrane is extended in one direction 
and crumpled in the other, with its crumpled bound- 
ary exhibiting statistics of a self-avoiding polymer. The 
usual crumpled-to-flat transition of an isotropic mem- 
brane is therefore generically split into two consecu- 
tive ones: crumpled-to-tubule and tubule-flat transitions. 
Fluctuations, nonlinear elasticity, and phase transitions 



into and out of the tubule phase have been studied 
extensivelv[lOj. with predictions dramatically confirmed 
by large-scale Monte Carlo simulations pi]. 

Motivated by unusual properties of nematic elastomers 
that have been recently discovered[H E fl E E E 
Hs[ . i.e. a spontaneous nematic order in an amorphous 
solid phase and resulting soft elasticity characterized by 
one vanishing shear modulus, we have previously stud- 
ied fluctuations and elasticity of a flat tethered mem- 
brane that exhibits spontaneous in-plane anisotropy, i.e. 
a nematic orderEl- Such spontaneous order guarantees 
the existence of a zero-energy in-plane deformation mode 
corresponding to a simultaneous reorientation of nematic 
and uniaxial distortion axes, and hence vanishing of an 
in-plane shear modulus. As in other "soft" fluctuating 
systems |2fJ, we found that elastic nonlinearities become 
qualitatively important in determining membrane's long- 
scale elasticity and conformations. However, due to the 
vanishing of an in-plane shear modulus, the elasticity 
of flat nematic elastomer membranes differ qualitatively 
from that of isotropic and that of explicitly anisotropic 
membranes: a renormalization analysis of the undulation 
nonlinearity shows that nematic-flat phase is character- 
ized by a conventional, length scale independent bending- 
rigidity modulus. E, HH 

While bulk nematic elastomers and gels have re- 
ceived considerable experimental and theoretical atten- 
tion, properties of equally interesting two-dimensional 
system, i.e. nematic elastomer membranes, have, to the 
best of our knowledge have not been explored in any de- 
tail. A few (rather idealized) experimental realizations 
of such spontaneously anisotropic nematic membranes 
can be envisioned. For example, one may prepare a ne- 
matic elastomer sheet in its isotropic phase and lower 
the temperature of the elastomer membrane into its ne- 



matic state. The chemical structure of liquid crystalline 
mesogenic units in such nematic elastomer should be cho- 
sen such that they prefer alignment parallel to the mem- 
brane surface. Liquid-crystalline polymers adsorbed onto 
a polymerized membrane can also develop a spontaneous 
nematic order at sufficiently high density 22], which in 
turn will act as a spontaneous elastic anisotropy to the 
membrane elasticity. The elastic properties of lipid bi- 
layers with spontaneous tilt order should be almost iden- 
tical to nematic elastomer membranes that we study in 
this paper |23| . Finally, other types of liquid crystalline 
orders may also be realized in polymerized membranes. 
For example, two dimensional smectic order may spon- 
taneously develops when temperature is further lowered 
in a nematically-ordered elastomer membrane. This can 
also be taken into account using a straightforward gen- 
eralization of our model. 

In this paper we continue to explore the physics of 
phantom |24f liquid-crystalline tethered membranes, fo- 
cusing on the overall phase diagram that results from the 
interplay of membrane's in-plane nematic and conforma- 
tional orders. As discussed in detail in Sec. IIHI we find 
that in the isotropic sector of the phase diagram, charac- 
terized by a vanishing in-plane nematic order, the system 
is identical to the usual isotropic tethered membranes, 
and therefore exhibits well studied isotropic-crumplcd 
(IC) and isotropic-flat (IF) phases. Development of ne- 
matic order generically leads to three additional phases: 
the nematic-crumpled (NC), nematic-tubule (NT) and 
nematic-flat (NF) phases. Different phases can be dis- 
tinguished by different expectation values of the nematic 
order parameter and the metric tensor. These results 
are summarized by a mean-field global phase-diagram in 

Fig.n 

We also study harmonic fluctuations within the most 
interesting nematic-tubule phase. As we discuss in 
Sec. IIVI and show schematically in Fig. [3 the geome- 
try and fluctuations of a nematic tubule, with an intrin- 
sic size LxL, are characterized by its average thickness 
Rg(L) and root-mean fluctuations h rms {L) about its en- 
ergetically preferred extended state, that scale as: 

R G {L)otL», hvms (L) oc L\ (1) 

where v and 9 are universal exponents. |ldj| For idealized 
"phantom" membranes v = 1/4 and — 1, but are ex- 
pected to be significantly modified by the self-avoidance 
and elastic nonlinearities. 10J These conformational prop- 
erties (or the absence thereof) can also be used as indica- 
tors of phase transitions into/out-of the nematic tubule 
phase, accompanied by usual thermodynamic singulari- 
ties, e.g., in the heat-capacity. For a more detailed dis- 
cussion of the scaling properties of polymerized tubules 
and transitions into/out-of the tubule phase, we refer the 
reader to Refllfl. 

We also derive a fully rotationally invariant nonlinear 
elasticity characterizing nematic tubules. We find that 
because of the spontaneous nature of the in-plane ne- 
matic order, the nonlinear elastic free energy of the ne- 
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FIG. 1: (Schematic mean-field phase diagram for a polymer- 
ized membrane with spontaneous in-plane nematic order. Pa- 
rameters r and j3 drive crumpling transition and isotropic- 
nematic transition respectively and are defined in SeclTTland 
Sec. IIHI Nematic-tubule phase is the stable phase in the 
shaded region. In the triangle region between three critical 
points Ci, C2 and C3, two phases are metastable. The dash 
line connecting Ci and C2 is a (schematic) first-order transi- 
tion line between isotropic/nematic-crumpled phases and the 
nematic-tubule phase. 
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FIG. 2: Conformational properties of the nematic tubule 
phase. Average thickness Rg and root-mean fluctuation forms 
are defined in Eqs. (I49H . Small rods, aligned on the mem- 
brane's surface schematically indicate existence of the in- 
plane nematic order. 



matic tubule differs qualitatively from the tubule phase of 
a permanently anisotropic membrane |10|. We leave the 
full analysis of tubule's elastic nonlinearities and ther- 
mal fluctuations to a future research. Here, we will have 
little to say on the very interesting and important ques- 
tion of the effects of heterogeneity, associated with the 
random cross-linking of elastomers. Based on a recent 
body of work on conventional liquid crystals confined in 
random environments |25, 26] and randomly polymerized 
membranes 27], we expect that at sufficiently long scales, 
even weak heterogeneities will lead to qualitatively im- 
portant modifications of some of the predictions made 
in this paper. Understanding these will be essential for 
a direct comparison with experiments!^]. We plan to 
explore these in the future. 

The body of this paper is organized as follows. In 
Section ^ we introduce a Landau model for a liquid- 



3 



crystalline tethered membrane, whose mean field phase 
diagram we map out in Section 11111 In Section II VI we 
focus on the most interesting phase, the nematic-tubule 
and study its thermal fluctuations at the harmonic level. 
In Section IIVBI we formally integrate out the nematic 
order parameter and derive an effective fully nonlinear 
and rotationally-invariant elastic free energy for the ne- 
matic tubule phase. 



II. MODEL OF NEMATIC ELASTOMER 
MEMBRANE 

Nematic elastomer membranes are direct generaliza- 
tion of permanently anisotropic tethered membranes 
studied in Ref. jlCj . with the in-plane elastic anisotropy 
chosen by a spontaneously, local nematic order. The 
latter plays the role of a new degree of freedom that 
fluctuates at a finite temperature. As for an ordinary 
tethered membrane, its geometric conformation (state) 
is described by a three-dimensional vector-field R(x), 
that specifies the position (embedding) of the reference- 
space mass point x in the three dimensional target space. 
Throughout this paper, we will use a lower-case bold-face 
characters, such as x = x a e a — xx + yy, to denote two- 
dimensional vectors in the reference space and upper-case 
bold-face character such as R = A" A" + YY + ZZ to de- 
note three-dimensional vectors in the target space, re- 
spectively spanned by a set of mutually orthogonal unit 
vectors, {ei = x, e-i = y} and {A, Y, Z}. 

The elastic free energy density of the membrane can be 
conveniently described as a function of a reference-space 
metric tensor g 

<9R <9R 

9ab ~ dx~ ' fe? {) 

inherited from the Euclidean metric Sij of the embed- 
ding space. By construction, g a b is explicitly a target- 
space scalar (i.e., rotationally invariant in the embed- 
ding space), rank-2 reference-space tensor, and is non- 
negative. The thermal expectation value of the metric 
tensor, (g) can be used to distinguish different phases of 
a tethered membrane. More specifically, two eigenvalues 
of (g) describe whether and how the membrane is ex- 
tended along two directions in its reference space. For 
example, (g) vanishes in the crumpled phase, and is pro- 
portional to the identity matrix in the flat phase, with the 
common eigenvalue describing the relative size of the flat 
membrane in the three dimensional target space. In the 
tubule phase of an intrinsically anisotropic membrane, 
(g) has only one positive eigenvalue, which corresponds 
to the fact that the membrane is extended along one di- 
rection and crumpled along others. Therefore (g) has all 
the properties of a physical order parameter. 

Nematic elastomer membranes are in addition charac- 
terized by a local in-plane nematic order parameter field 



Q(x). Characterizing membrane's spontaneous elastic 
anisotropy, Q describes an intrinsic property of the mem- 
brane and is a second-rank symmetric traceless tensor in 
the two dimensional reference space, with components 
Qab, Q, b = x,y and ^ Q Q aa = 0, and is a target-space 
scalar 29j. Landau free energy density for our model of a 
nematic elastomer membrane is a rotationally-invariant 
(both in the reference and target spaces) function of g, 

£: 

fig, £] = fe [g] + fn [£] + fen [g, Q], (3) 

where f e is the elastic free energy density for an isotropic 
polymerized membrane 0, 

U [g] = \ (V 2 R) 2 + r'Trg + \ (Trg) 2 + M Tr £ 2 . (4) 

/„ is the Landau-deGennes free energy density for a 
nematic-isotropic transition in two dimensions 

/n = f (V£) 2 + fTr£ 2 + |(Tr£ 2 ) 2 . (5) 

It reflects the special feature of nematic-isotropic transi- 
tion in two dimensions, namely, the absence of the usual 
cubic invariant TrQ 3 that vanishes identically. This al- 
lows nematic-isotropic transitions in two dimensions to 
be continuous [30j. 

The last part of /, Eq. is the coupling between ne- 
matic order and the elastic degrees of freedom, which, to 
lowest order, is given by: 

f en = -aTrQg = -o/Lr Qg, (6) 

where g = g— {Trg) 7/2 is the traceless part of the metric 
tensor. In the presence of such tensorial nemato-elastic 
coupling, nematic order Q acts as a uniaxial shear stress 
and therefore renders the in-plane elasticity anisotropic 
inside a nematic phase. Concomitantly, an elastic distor- 
tion characterized by an anisotropic metric tensor acts 
as a quadrupolar field inducing nematic order. As usual 
in Landau description, although other higher order cou- 
plings exist, such as for example "/TrgTrQ 2 , they do not 
lead to any new qualitative effects, and, furthermore, 
for weak nematic order, are quantitatively subdominant. 
Their primary effect is to slightly change the locus of 
phase boundaries without changing their topology. We 
will therefore ignore them in this paper. 

We can also formally and perturbatively integrate out 
the nematic order parameter Q in Eq. @ and obtain an 

effective elastic free energy density formulated purely in 
terms of the metric tensor. Up to fourth order in <?, it is 
given by 

fes[g] = ^V 2 r) 2 + im 1 + M (Tr| 2 )+^Trg) 2 

+ B(Tr| 2 ) 2 + CTr 5 Tr| 2 , (7) 
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and agrees with the purely elastic model used in our pre- 
vious study of the nematic-flat phase 19j. Although such 
effective model can also be used to study the phase dia- 
gram and other phases, it cannot distinguish the nematic- 
crumpled and the isotropic-crumpled phases, since both 
are characterized by a vanishing metric tensor. There- 
fore, depending on the equilibrium value of the metric 
tensor, controlled by model parameters, / e ff predicts a 
crumpled phase, a spontaneous tubule phase (nematic- 
tubule phase), an isotropic-flat phase, and a spontaneous 
anisotropic flat phase (nematic-flat phase). Here, the 
usual crumpled-to-flat transition is determined by chang- 
ing the sign of t, while the crumpled-to-tubule transition 
is driven by a sign of u oc /i. This purely elastic approach 
has the advantage of simplicity, stemming from absence 
of the nematic order parameter. In Sec II VI we will use 
this effective elastic free energy to study the nonlinear 
elasticity in the nematic-tubule phase. 



III. MEAN-FIELD PHASE DIAGRAM 

An isotropic polymerized membrane with elastic free 
energy density, Eq. (gj undergoes a continuous^ crum- 
pling transition as r' is tuned from positive to negative Q. 
On the other hand, /„, Eq. exhibits a second-order 
isotropic- nematic transition as (3' changes sign[3(|. When 
both r' and j3' are small, a competition between these two 
types of order becomes important. 

To work out the mean-field phase diagram, we mini- 
mize the free energy density Eq. (0) over the metric tensor 
g and the nematic order parameter Q. To this end, it is 
convenient to rescale both order parameters g and Q as 
well as the total free energy Eq. |21 as following, 



where 



a 



a 2 



Q. 



16\/ fi 3 v — 

/• 



256/i 2 v 



(8a) 
(8b) 
(8c) 



Furthermore, it is clear that the free energy is minimized 
by uniform nematic order and metric tensor. Therefore 
we can ignore the two higher derivative terms k, K ap- 
pearing in Eqs. Q and JSJ, respectively. The resulting 
rescaled total free energy density has the following simple 
form: 



f\9,Q] 



(Tr(g +rl- Q)f + Tr(g + rl- Qf 



+ 7tW + /5) 2 



(9) 







a 2 (A + /i)' 
1 

a 2 ' 
A 



(10a) 
(10b) 
(10c) 



and J is the 2x2 identity matrix. Stability of the elastic 
free energy Eq. |@} requires fi > and > 0, therefore 
e>-l. 

The linear coupling Eq. 10 between g and Q prefers 

same "orientation" (alignment of their eigenvectors) for 
these two tensors. Therefore in the ground state, where 
g and Q minimize the total free energy, Eq. @, these 
two tensors commute. Using the coordinate system in 
which g and Q are simultaneously diagonal, they can be 

parameterized as: 



9i 
g 2 



0=1 s 

M \ -S 



(11) 



where g\ and gi are two non-negative eigenvalues of the 
metric tensor g, and S a non- negative magnitude of the 
nematic order. Substituting these into free energy /, 
Eq. we find that it can be written as 



/ = jj/ff [01.02, 5] + ^fs[S], 



(12) 



where 



f g [gi,92,S] = (2+ e ) (( 9l - ai ) 2 + (g 2 - a 2 ) 2 ) 



+ 2e(gi-ai)(g 2 -a 2 ), 

= (1 + e)(gi + g 2 - ai - a 2 ) 2 , 

+ (ffl - 92 - Oil + «2) 2 , 

fs[S) = (S 2 +p) 2 , 



and 



ai = (S-r), 
a 2 = -(S + r), 



(13a) 

(13b) 
(13c) 



(14a) 
(14b) 



with ai > a>2, and f g , fs non- negative definite. 

To obtain a mean-field phase diagram, we need to min- 
imize /, Eq. 1 121 over nonnegative gx, g 2 and S for given r 
and p. Before doing that, however, it is useful to classify 
all possible phases according to ground state values of 
<7i, g 2 and S. It is easy to see that for vanishing nematic 
order, the free energy density, Eq. © reduces to that 
of an isotropic polymerized membrane, i.e. Eq. Q. An 
isotropic polymerized membrane is either in the isotropic- 
crumpled phase where gi = g 2 = 0, or in the isotropic-flat 
phase, where gi — g 2 > 0. In other words, a state with 
9i 92 and S — cannot be a stable ground state. Con- 
versely, anisotropy in the metric tensor, g\ ^ g 2 acts as 
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an external field on the nematic order parameter Q, in- 
evitably inducing finite nematic order. Therefore S must 
be finite whenever g\ ^ g 2 . Similarly, a nonzero nematic 
order acts as an external shear stress tensor, which always 
induces a finite shear strain in a flat membrane (Hook's 
law), leading to elastic anisotropy with gi ^ g 2 - There- 
fore g± and g 2 cannot be identical in the presence of a 
finite nematic order, S. That is, if S ^ 0, the membrane 
can either be in a nematic-crumpled phase (g± = g 2 = 0) , 
or a nematic-tubule phase (gi > g 2 = 0), or a nematic- 
flat phase (31 > g 2 ^ 0). We list the resulting five mem- 
brane phases in Table. [I] 



gi, gi and S 


S = 




gi = gi = 


Isotropic Crumpled 


Nematic Crumpled 


gi = gi > 


Isotropic Flat 


Unstable 


gi > gi > 


Unstable 


Nematic Flat 


gi > gi = 


Unstable 


Nematic Tubule 



TABLE I: Possible phases of a nematic elastomer membrane. 

We now minimize the rescaled free energy, Eq. i|12|) 
over gi, g 2 and S. It is easy to see that if (3 > 0, fs[S] 
in Eq. (jl3c|l is minimized by 



5 = 0. 



(15) 



Substituting this into Eq. (|13af> . we find that 
f g [di, a2, S = 0] is minimized by 



,9i = .92 



0, for r > 
/•|, for r < 0. 



(16) 



which, as discussed above corresponds to the isotropic- 
crumpled and the isotropic-flat phase, respectively. This 
situation is illustrated in the (3 — r plane in Fig. ^ 

For (3 < 0, our strategy will be to determine the region 
of (meta-)stability for various phases. We first minimize 
fg[gii9i,S\, defined by Eq. (|13af> . for given a\ and a 2 , 
i.e. for given S and r. It is clear that for fixed values 
of ot\ and a 2l a contour f g = Const, is an ellipse in the 
(51 j 92) plane centered at (ai,a 2 ). Its principle axes are 
along two diagonal (1, ±1) directions, as shown in Fig. [3] 
Depending on the values of ol\ and a 2 , i.e. the center of 
the ellipse, the minimizer of f g [gi, g 2 , S] falls into three 
categories: 

1. ai > a 2 > : 

The ellipses are centered within the first quadrant. 
It is easy to see that f g is minimized by 



(9*i^9* 2 ) = {ui,a 2 ) = (S-r,-S-r), 



(17) 



with a minimal value zero. This corresponds to a 
nematic-flat phase. Substituting these g\ and g 2 
into Eq. O, we obtain: 



(18) 





g 2 ^~^<gi= £2 
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FIG. 3: Contour plot of f g — const., an ellipse centered at 
(0:1,0:2) in the ((71,32) plane, where ai < 0:2, with principle 
axes along diagonal (1, 1) direction. As discussed in the text, 
in case 1, f g is minimized by the center of the ellipse, that cor- 
responds to a nematic-flat phase. In case 2, it is minimized by 
the point T, where the ellipse is tangent with the positive gi 
axis. This corresponds to a nematic-tubule phase. In case 3, 
f g is minimized by the origin g\ — gi = 0, which corresponds 
to a (nematic or isotropic) crumpled phase. 



which, when minimized over S, gives the magni- 
tude of nematic order parameter in the nematic-flat 
phase: 



s = V\P\- 



(19) 



The condition of non- negativity of g\^ 2 and g 2 = a 2 , 
gives the condition of stability of the nematic-flat 
phase 



9*2 = ~(S + r) 

= -(r+Vi)>0, 



(20) 
(21) 



which leads to the phase boundary for the nematic- 
flat phase illustrated in Fig. ^ 



2. ai 



2+£ 



a 2 > and a 2 < 0: 



In this case, the centers of elliptic free energy con- 
tours are located outside the first quadrant, and 
therefore do not correspond to physical solution of 
the metric tensor, characterized by gi >2 > 0. Inde- 
pendent of ai, however, there is always one elliptic 
contour, which is tangent to the positive g\ axis as 
illustrated in Fig. 0. It is straightforward to see 
that this tangent point, (gl > 0,g 2 = 0) is the min- 
imum of f g in the physical space of metrics, which 
corresponds to the nematic-tubule phase. At this 
tangent point, the partial derivative of f g with re- 
spect to (71 vanishes, i.e. 



dgi 



(95%0) 



2(2 + e)Gtf-ai 



ea 2 



) = 0, (22) 
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which leads to 



9l 



ai 



ea 2 



2 + e 
2(S-(l + e)r) 
2 + e 



(23) 



Substituting this g\ and g% = into Eq. ljT2l . the 
rescaled total elastic free energy reduces to 

The equilibrium value of 5* is determined by mini- 
mizing this / with constraints 



a 2 < 
> 



5 > -r, 

5 > (1 + e)r. 



(25a) 
(25b) 



For a negative r, it is easy to see that as long as 
r < \f\j3\, there is always a positive value of S that 
both minimizes Eq. (|24|l and satisfies constraints, 
Eqs. - The case 01 positive r is more subtle. By 
a technically involved, but conceptually straightfor- 
ward calculation, we find that for (3 < — 2— (i.e. 
to the left of point C\ in Fig.QJ, a nematic tubule 
state is stable as long as (1 + e)r > \/|2 + /3|. For 
-2(1 + e)/(2 + e) > /3 > -2 - (i.e. to the right 
of Ci, but to the left of C 2 in Fig. QJ, we find that 
a metastable nematic-tubule state exists as long as 



2(1 



< 



(26) 



And finally for -2(1 + e)/(2 + e) < (3 < (i.e. to 
the right of C2 in Fig. [TJ, nematic tubule phase is 
stable if and only if r < 0. 

a>i + 2^7 a 2 < and a 2 < 0: 

In this case, there is no ellipse tangent with the pos- 
itive gi axis, i.e. the nematic-tubule phase is not 
stable. There is however always one ellipse pass- 
ing through the origin, g* = g\ = (see Fig. |3J) 
that minimizes f g . Depending on the value of S, 
the corresponding phase may be nematic-crumpled 
or isotropic-crumpled. Setting g\ = g 2 = in 
Eq. flUSjl. we obtain 



/ = Is- 

1 2 



+ (2 + (3)S 2 + ^/3 2 + 2(1 + e)r 2 



(27) 



where we have used Eq. I|14afl and Eq. (|14b[) . Min- 
imizing it over S, we find the equilibrium value of 
nematic order parameter 



S = 



3 



if (3 + 2 > 
if (3 + 2 < 



(28) 



Considering the self-consistent condition that de- 
fines case 3, as well as the definition of a\ and a 2 , 



we find that for negative /3, the isotropic-crumpled 
phase is stable if and only if 



/3 + 2 > and r > 0. 



(29) 



On the other hand the nematic-crumpled phase is 
stable if and only if 



(3 + 2 < and < a/1/3 + 2| < (1 + e)r. (30) 

We have thereby determined the region of meta- 
stability for all five phases admitted by the nematic- 
elastomer membrane. These results are summarized in 
Fig-H I n the ((3, r) plane the small triangular region be- 
tween three critical points C\ , C 2 and C3 is particularly 
interesting. This region is divided into two parts by the 
vertical line (3 = —2 passing through C3. In the left part, 
both nematic-tubule and nematic-crumpled phases are 
metastable. In the right part, both nematic-tubule and 
isotropic-crumpled phase are metastable. Consequently 
there exists a first-order transition line between C\ and 
C 2 . This line is illustrated schematically by a dashed 
curve in Fig. ^ It connects a second-order nematic- 
tubule — to — nematic-crumpled transition boundary on 
the left to a second-order nematic-tubule — to — isotropic- 
crumpled phase boundary on the right. The exact posi- 
tion of this first-order line, which can in principle be de- 
termined from comparing free energy of different phases, 
is not essential to our discussion. Outside the triangular 
region, all phase transitions are continuous in our mean- 
field analysis. 

For completeness, in Table. ITTlwe list the coordinates of 
the three critical points, C\, C 2 and C3 in the ((3, r) plane. 
It is interesting to note that distances between any two of 
these points, and therefore the triangular area, vanishes 
as e becomes large, corresponding to an incompressible 
limit of vanishing shear-to-bulk moduli ratio. Due to the 
same reason, in the limit of small shear modulus, the 
discontinuity of first-order transition line (jump of order 
parameter or other physical quantities) vanishes. 





Ci 


c 2 


c 3 








-2 


r 


1 

(l + e)V2 + f 









TABLE II: Coordinates of three critical points Ci, C2 and 
Cs. 



IV. PHANTOM NEMATIC-TUBULE PHASE 

Among all phases that we have identified in the pro- 
ceeding section, the nematic tubule phase is the most 
interesting one. As we have just found, it is character- 
ized by a nonzero nematic order and a metric tensor with 
one positive eigenvalue, i.e. 



£0 = 





-s 



9i 




(31) 
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Even though the ground state configuration of this ne- 
matic tubule phase is identical to that of the tubule 
phase studied in Ref . , as we will show in this section, 
its elasticity and thermal fluctuations differ qualitatively, 
due to the fact that rotational symmetry is spontaneously 
broken, both in the reference space and in the target 
space. 

A. Harmonic Fluctuations 

We first study harmonic thermal fluctuations in the 
nematic tubule phase. Without loss of generality we pick 
the intrinsic (reference space) membrane axes such that 
in the ground state the x-axis coincide with the extended 
direction of tubule, with tubule therefore crumpled along 
the y-axis. Similarly, our choice of the target space co- 
ordinate system orients the tubule along X-axis of the 
embedding space. With this choice, the tubule ground- 
state conformation is given by 

R (r) = QxX, (32) 

where g\ — C, 2 can be calculated by minimizing the total 
free energy Eq. (JSJ. 

To study thermal fluctuations about the nematic- 
tubule ground state, we follow the parameterization in- 
troduced in Ref. [l(| for the study of an intrinsically- 
anisotropic tubule. A deviation away from the ground 
state can be parameterized by a one-dimensional (scalar) 
phonon field u(r) and a two-dimension undulation field 
h(x), defined by 

R(r) = (C.x + u(r))A + h(r). (33) 

The corresponding metric tensor, Eq. is then given 
by 

g° xx = (C+M' + ^h) 2 , (34a) 
g% = (C + d x u)(d y u) + (d x h) ■ (d y h), (34b) 

9°yy = (d y uf + (d y hf. (34C) 

The Lagrange strain tensor is defined as half the devia- 
tion of the metric tensor from its ground state value g , 
i.e. 

Uab=^(9ab- 9ab)- (35) 

Its components are given by 

u xx = C (d x u) + ^(d x uf + ^{d x h) 2 , (36a) 
u xy = -(C + d x u)(d y u) + -(d x h) ■ (dyh), (36b) 

«i* = + (36c) 

Fluctuations of the nematic order parameter Q about 
its ground-state value Q can be parameterized by two 



variables r and a defined by: 

£ = £ + 5 £=( 5 r -s-r)- w 

They correspond to longitudinal and transverse excita- 
tions of the nematic order, respectively. Their physical 
significance can be readily seen by comparing Eq. I|37|l 
with the more usual parameterization of nematic fluctu- 
ation in terms of rotation angle 9 and magnitude fluctu- 
ation SS: 

2=<^>(s 2 2 r s ii fl > < 38 > 

We find that up to linear order of SS and 8, 

a « 2S0, t w SS. (39) 

Therefore the longitudinal component r describes fluctu- 
ations in the magnitude of nematic anisotropy about S, 
while the transverse component a is proportional to the 
rotation of the nematic director. 

When expressed in terms of these scalar fields and the 
Lagrange strain tensor, the nematic-tubule free energy, 
Eq. becomes 

/[£,£] = ^-^{2r + C 2 + 2u xx + 2u yy ) 2 , 

+ ^[f3 + 2a 2 + 2(5 + t) 2 ) 2 + 2(<t -2u xy f 
+ 2 (S - ( 2 /2 + T -u xx + uyyf . (40) 

The condition that g and Q , as given by Eq. 131fl . 
minimize the elastic free energy density f[g,Q] guaran- 
tees that in Eq. 1)40(1 . terms that are linear in r or u xx 
strictly vanish, and thereby leads to two equations that 
determine the ground-state order parameters 5* and C : 

~2S + 2r (1 + e) + (2 + e) ( 2 = 0, (41a) 
4S 3 + 2S {l+{3)-( 2 = 0. (41b) 

Finally, eliminating the strain components in favor 
of the phonon (u) and height undulation (h) fields using 
Eqs. (|36|l . and keeping only upto quadratic terms, we 
obtain the harmonic free energy density of the nematic- 
tubule: 

f[a,r,u,h] = /[£,£] -/[£ ,£J 

= K a {Va) 2 + C a (a- a<J dyu) 2 

+ K t {\7t) 2 +C t {t- a T d x u) 2 

+ B u (d x u) 2 + K u (d y 2 u) 2 

+ B h (d y h) 2 + K h (d x 2 h) 2 1 (42) 



8 



where 



C T 

B u 
B h 



2 (1 + 2S 2 
2 (1 + 6S 2 

c 



■0), 



1 



2S 2 

c 



■0' 



1 + 6S 2 + 13' 

2 (2 + e) C 2 , 
4S-2( 2 . 



(43a) 
(43b) 

(43c) 

(43d) 

(43e) 
(43f) 



There are several salient features of the nematic-tubule 
free energy density, Eq. H42JI . that are worth pointing out. 
Firstly, terms that are first-derivative of h with respect 
to the ordered direction, x, namely (d x h) 2 , do not appear 
in /. Small d x h corresponds to an infinitesimal rotation 
in the embedding space of the tubule's extension axis, 
that we have arbitrarily chosen to point along X and 
therefore must not cost any elastic free energy. There- 
fore absence of the (d x h) 2 term is a direct consequence 
of the underlying rotational invariance of the embedding 
space, that a tubule, extended along X, spontaneously 
breaks. Secondly, / exhibits a linear coupling of the 
longitudinal fluctuation of the nematic-order parameter, 
r to the longitudinal strain, d x u. This corresponds to the 
linear coupling between Q and the metric tensor g in the 

original model Eq. © and encodes the fact that nematic 
order induces uniaxial strain and vice versa. Lastly and 
most importantly, the transverse fluctuation of the ne- 
matic order, a is "minimally" |32j coupled to the trans- 
verse strain component d y u, with the coefficient of the 
linear coupling o~d y u precisely such that it forms a com- 
plete square with two other terms a 2 and (d y u) 2 . This 
ensures that a small and nonzero transverse strain d y u 
(corresponding to a reorientation of the uniaxial strain 
axis) can be completely compensated by a uniform per- 
turbation a, which corresponds to a global rotation of the 
nematic director within the reference space, so that the 
overall free energy change is exactly zero. This property 
is a consequence of the underlying rotational invariance 
of the membrane's reference space, spontaneously broken 
by the nematic order, and is therefore expected on gen- 
eral symmetry principles. One consequence of such cou- 
pling is that integration over the nematic director field, tr, 
leads to an effective elastic theory for the nematic-tubule, 
in which the shear term (d y u) 2 is absent. The resulting 
extra "softness" is the essential feature that qualitatively 
distinguishes a nematic-tubule studied here from a tubule 
in an explicitly anisotropic polymerized membrane stud- 
ied by Radzihovsky and Toner. 0] It reflects the fact 
that in a nematic-tubule the rotational symmetry in the 
reference space is only spontaneously broken, and there- 
fore a simultaneous rotation of the nematic order and the 
metric tensor in the reference space does not cost energy. 
For mechanical stability we have therefore added to the 
free energy Eq. (|42|l two curvature terms, K u and Kh, 
that would otherwise be unimportant for smooth defor- 



mations. 

To study thermal undulations of the nematic-tubule it 
is convenient to integrate out fluctuations of nematic or- 
der parameter, r and a, in a direct analogy of integrating 
out the nematic director to study properties of the smec- 
tic phase. This leads to an effective elastic Hamiltonian 
for the nematic-tubule (in harmonic approximation) 

H eS [u,h] = ~ J dxdy[B u {d x u) 2 +K u (8 y 2 u) 2 

+ B h (d v h) 2 + K h (d x 2 h) 2 ] , (44a) 

= \ J dq x dq y [(B u q 2 x +K u q 4 y )\u(q)\ 2 

+ (B h ql + K hq i)\h(q)\ 2 ] : (44b) 

where in the last equation we have also expressed the 
elastic Hamiltonian in terms of the Fourier components 
defined by: 



«(q) = Jdxdye-^ „(r), 
h(q) = / dxdye- tcl r h(r). 



(45a) 
(45b) 



From -ffgff we observe that at the harmonic level, fluc- 
tuations of h field and the phonon field u are decoupled 
from each other. Furthermore, the in-plane phonon-field 
u elasticity is identical to that of extensively studied two- 
dimensional smectic liquid crystal with layers extended 
along Therefore we can immediately adopt 

known results for a two-dimensional smectic to fluctu- 
ations of the nematic-tubule phase. |l9i l33j For example, 
for a membrane with intrinsic size L x x L y , the bulk 
contribution to mean-squared fluctuations (see below for 
contribution from other modes) for u is given by: 



(u 2 ) 



1 



i ) w ^-^ipu{2Tra u L x /L 2 ) 
B„ V Zira,, 



(46) 



1 I L y 
(2tt) 2 B u a u ' 



if Ly > 2na u L x ; 
if L 2 < 2na u L x 



Therefore, thermal fluctuations of the phonon field, u 
grow without bound with membrane's intrinsic size (L x , 
L y ). This growth of thermal fluctuations of the u field is 
a signature of the qualitative importance of the hitherto 
neglected nonlinear elastic terms. 

From Eq. I|42|l we observe that at long length scales 
(beyond \J K a /C a ) the nematic orientation field, a is 
locked by the "minimal coupling" term to transverse 
strain fluctuations, d y u. This, together with dimensional 
analysis allows us to relate the asymptotic scaling of the 
nematic director field, a to that of the phonon field u 
through 



cx ^ (u 2 



(47) 
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Despite a strong growth of (it 2 ) appearing in Eq. (|46|l . 
we therefore find that the in-plane nematic orientational 
fluctuations remains finite in the thermodynamic limit. 

From -ffgff we observe that the undulation field, h, 
is also governed by the same two-dimensional smectic 
anisotropic elasticity, except that the roles of x and y axes 
are interchanged. At the harmonic level, it is identical to 
that of a permanently anisotropic membrane studied in 
Ref. As discussed in great detail there, two related 
quantities can be defined to describe the tubule geome- 
try at a finite temperature, tubule radius of gyration, Rq 
(thickness) and tubule roughness, h rms as follows: 

R 2 G = (Ih^-hOr,^)! 2 ), (48a) 
h 2 rms = (\h(0,y)~h(L x ,y)\ 2 ). (48b) 

As indicated in Fig [21 Rq measures the radius of a typ- 
ical cross-section of the tubule that is perpendicular to 
the nematic order and tubule axis, while h rms measures 
wandering of the tubule transverse its backbone. As 
discussed by Radzihovsky and Toner 10J fluctuations of 
these two quantities are strongly affected by the so-called 
zero modes, that correspond to Fourier modes of the un- 
dulation field h with one vanishing wave vector compo- 
nent. Adapting these results to our case, we find: 

Rg{L x , L y ) = L^ 4 Sr(L x / \/Ly), (49a) 
h rms (L x ,L y ) = L 1 J 2 S h {L x l^fT y ), (49b) 

where two crossover functions Sr and Sh satisfy 

S R (x) - X ^°> (50) 

I constant x — > oo; 

„ , . I constant x — > 0, ._.„, 
Sh{x) -> < (51) 
I a; x — + oo. 

For an approximately square size membrane (L x ~ L x ~ 
L), we find that h rms scales linearly with L, while Rq cx 
L 1 / 4 . Finally, we note that quantities similar to Eq. Ij48(l 
can also be defined for the fluctuations of the phonon field 
u. Their asymptotic forms are also given by Eq. (|49|l . 
with the role of x and y axes interchanged. 

As was found for explicitly-anisotropic tubules(lfT|. we 
expect that inclusion of the so-far ignored self-avoidance 
and elastic nonlincarity (both from in-plane phonon as 
well as height undulation) will swell the nematic-tubule 
and further stabilize it against thermal fluctuations. In 
order to perform such a study, we will need a formalism 
of elastic free energy for the nematic tubule that includes 
all relevant nonlinearities. We derive this nonlinear elas- 
ticity in the next subsection and leave the study of it to 
a future research. 



B. Nonlinear theory for nematic-tubule 

Our goal now is to derive the fully rotationally- 
invariant (in both the reference and target spaces) free 



energy describing nematic-tubule. This turns out to be 
simpler within the purely elastic formulation, in which, as 
discussed in Sec|n] the in-plane nematic orientational de- 
grees of freedom have been formally integrated out. The 
analysis here, then parallels our work on nonlinear elas- 
ticity of three-dimensional, bulk nematic elastomers. |l7j 
In such purely elastic approach nematic-tubule's free 
energy is a functional of the deviation of metric tensor 
g from its equilibrium values go- Since the elastic free 
energy density should be invariant under arbitrary rota- 
tions in both the reference and the target spaces, for ho- 
mogeneous deformations and in two dimensions it must 
be a function of only two scalar functions of the metric 
tensor j34|: 

51 = Trg, (52a) 

5 2 = Trg 2 . (52b) 

In the nematic-tubule phase, we can expand the effec- 
tive elastic free energy density around the ground state 
corresponding to the metric tensor go , Eq. I|31|) : 

f eS [Si,S 2 ] = f eS [Sl,S%]+A 1 5S 1 + A 2 5S 2 

+ ^BxSS 2 + Bi2<5Si5S 2 + \b 2 8S\ 
+ 0(SSf), (53) 

where 

551 = Trg - Trg = 2 Tr& (54a) 

55 2 = Trg 2 - Tig 2 = 4 Trg u + 4Tr^ 2 . (54b) 

As can be checked a posteriori, terms of order 
SSf and higher are irrelevant at long-scales in the 
renormalization-group sense and can therefore be omit- 
ted, as we have done above. 

Using Eq. (|54|l . we may further express the elastic free 
energy density in terms of the Lagrange strain tensor. To 
proceed, however, we note the following relations between 
SSi and 5S 2 : 

SS1SS2 = 2( 2 5S 2 + irrelevant 

= 8C 2 (Trj|) 2 + irrelevant, (55a) 

SSi = 4 C 4 6S 2 + irrelevant 

= 16C 4 (Tr^) 2 + irrelevant. (55b) 

Substituting these relations, together with Eqs. (|54|l into 
Eq. H53|) . we obtain the following nonlinear elastic free 
energy density for the nematic tubule: 

fM = f C ff[Si,s 2 }- Msls"} 

= 2{A X + 2( 2 A 2 )Trg + 4A 2 {-( 2 u yy + Tr^ 2 ) 
+ 2(S 1 +4C 2 S 12 +4C 4 S 2 )(Tr i ) 2 . (56) 

Using the fact that g^, Eq. l(3*T|) is a true ground state, 
together with the observation, from Eqs. (|36fl . that Tru 
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and u yy are, respectively, linear and quadratic in the dis- 
placement fields, u and h, we find that the coefficient of 
Tru identically vanishes, namely 

Ai + 2( 2 A 2 = 0. (57) 

To express the elastic free energy in terms of the phonon 
field u and the undulation field h, we furthermore note 
the following two identities: 

Trj| = C^M+i(Vii) 2 + i(Vh) 2 , (58a) 

-( 2 u yy + Tm 2 = (Ca^+i(Vw) 2 + i(Vh) 2 ) 2 

- i (Cd y h + (d x u)(d y h) - (dyu)(d x h)) 2 

+ irrelevant (58b) 

Finally, in order to further simplify the notation, we 
rescale the coordinates in the reference space according 
to: 

x^C 1 ^, (59) 

thereby obtaining our final expression for the nonlinear 
elastic free energy density of a nematic tubule: 

f[u,h] = l -B u [{d x u)+ l -{Vuf + i(Vh) 2 ] 2 

+ ^B h [dyh(l + d xU ) ~ (d v u)(d x h)] 2 

+ ^K u (d v 2 u) 2 + ^K h (d x 2 h) 2 , (60) 
with the effective bulk elastic moduli given by: 

B u = 8C 4 A 2 +4C 4 (Bi+4C 2 Bi2+4C 4 B 2 ), (61a) 
B h = -4C 4 A 2 . (61b) 

Despite our ignoring of a number of (at long-scales) 
subdominant terms, it is easy to verify that our final 
expression for the elastic free energy density, Eq. (|60|l . 
is actually exactly invariant under arbitrary rotations in 
both the reference and the target spaces. To see this 
we first note that tubule's reference ground-state confor- 
mation, using the rescaled coordinate system Ea. (|59|l . is 
given by 

R (x)=Xx. (62) 

Under arbitrary rotation both in the reference space and 
in the target space, it becomes 

R(x) = N {h ■ x) = X(x + uo(x)) + h (x), (63) 

characterized by target and reference-space unit vectors, 
Nq a n d ho, respectively. This situation is illustrated in 
Fig. (|IV B|l . Without loss of generality, we may parame- 
terize these two unit vectors as 



y 




FIG. 4: A ground state conformation of a nematic tubule can 
be parameterized by two unit vectors no and iVo, Eq. 163L 
that respectively describe the orientation of the tubule in the 
reference and target spaces. In the reference state Eq. 162L 
ho and No are chosen to be x and X respectively. 



From Eq. I)63[) . it is easy to see that relative to the refer- 
ence state Ro, Eg. 1)62)) . this pure rotation corresponds to 
a phonon field ito( x ) and an undulation field ho(x) given 
by 

ito(x) = £a;(cos 8 cos (j> — 1) + Qy sin cos 0, (65a) 
ho(x) = sm<f)(xcos6 + ysm8). (65b) 

Substituting these two equations into f[u, h], EqEOI we 
find, that, indeed, as demanded by rotational invariance, 
the nonlinear elastic free energy strictly vanishes for an 
arbitrary rotation in both the reference and target spaces. 



V. CONCLUSION 

In this paper, we have presented a Landau theory for 
polymerized membrane with spontaneous in-plane ne- 
matic order. We have analyzed the mean-field phase di- 
agram and studied harmonic thermal fluctuations in the 
nematic-tubule phase. We have also derived a nonlinear 
elastic free energy for the nematic-tubule, which will be 
an essential starting point for future analysis of elastic 
nonlinearities and self-avoiding interaction that we have 
so far ignored. Such a study is necessary in order to un- 
derstand the physics of any realistic membrane exhibiting 
spontaneous nematic tubule phase. 



Acknowledgments 



A^o = X cos <f> + Y sin 0, (64a) 
h = icos^ + ysin6*. (64b) 



We thank T. C. Lubensky for discussion and acknowl- 
edge support by the NSF MRSEC DMR-0213918 (LR), 



11 



DMR-0321848 (LR), David and Lucile Packard Founda- partment of Physics and KITP at UC Santa Barbara, 
tion (LR), as well as NSF DMR02-05858 (XX)and DOE where part of this work was done, for hospitality. 
DEFG02-91ER45439(XX). We also thank Harvard De- 



[1] Also referred as fixed connectivity, tethered, or crystal 
membrane in various references. 

[2] For a review, see Statistical Mechanics of Membranes and 
Interfaces, edited by D. R. Nelson, T. Piran, and S. Wein- 
berg (World Scientific, Singapore, 1989). 

[3] D. R. Nelson and L. Peliti, J. Phys. (Pans) 48, 
1085(1987). 

[4] J. A. Aronovitz and T. C. Lubensky, Phys. Rev. Lett. 
60, 2634(1988); J. A. Aronovitz, L. Golubovic, and 
T. C. Lubensky, J. Phys. (Pans) 50, 609(1989). 

[5] P. LeDoussal and L. Radzihovsky, Phys. Rev. Lett. 69, 
1209(1992). 

[6] Y. Kantor and D. R. Nelson, Phys. Rev. Lett. 58, 
2774(1987). 

[7] M. Paczuski, M. Kardar, and D. R. Nelson, Phys. Rev. 

Lett. 60, 2638(1988). 
[8] E. Quitter, F. David, S. Leibler, and L. Peliti, Phys. Rev. 

Lett. 61, 2949(1988). 
[9] L. Radzihovsky, to appear in the 2nd Edition of the Sta- 
tistical Mechanics of Membranes and Surface, edited by 
D. R. Nelson, T. Piran, and S. Weinberg, World Scien- 
tific, Singapore. 
[10] L. Radzihovsky and J. Toner, Phys. Rev. Lett. 75, 

4752(1995); Phys. Rev. E 57, 1832 (1998). 
[11] M. Bowick, M. Falcioni, and G. Thorleifsson, Phys. Rev. 

Lett. 79, 885(1997); cond-mat/9705059 
[12] L. Golubovic and T. C. Lubensky, Phys. Rev. Lett. 63, 
1082(1989). 

[13] M. Warner and E. M. Terentjev, Prog. Polym. Sci. 21, 

853(1996), and the reference therein. 
[14] E. M. Terenjev, J. Phys. Cond. Mat. 11, R239(1999). 
[15] H. Finkelmann, I. Kundler, E.M. Terentjev, and M. 

Warner, J. Phys. II 7, 1059(1997); GC. Verwey, M. 

Warner, and E.M. Terentjev, J. Phys. II (France) 6, 

1273-1290(1996); M. Warner, J. Mech. Phys. solids 47, 

1355(1999). 

[16] T. C. Lubensky, R. Mukhopadhyay, L. Radzihovsky, 
and X. Xing, P hys. Rev. E 66, 011702(2002); 
cond-mat/0112095 

[17] X. Xing and L. Radzihovsky, Europhys. Lett. 61, 769 
(2003); Phys. Rev. Lett. 90, 168301(2003). 

[18] O. Stenull and T. C. Lubensky, Europhys. Lett. 61, 776 
(2003); Phys. Rev. E 69, 021807(2004). 

[19] X. Xing, R. Mukhopadhyay, T. C. Lubensky, and 
L. Radzihovsky, Phys. Rev. E 68, 021108 (2003). 

[20] Similar phenomena take place in other "soft" anisotropic 
systems, such as smectic (25| and columnar liquid crys- 
tal phases |25l. vortex lattices in putative magnetic 
^sup erconductors|26 | . and bulk nematic elastomers|l6l.ll^. 

[21] Stability of an isotropic-flat phase at finite temperature is 
ensured by length scale dependent bending modulus that 
is infinitely enhanced by thermal fluctuations. Therefore 
this finding of a finite, length-scale independent bending- 
rigidity modulus for flat nematic elastomer membranes 
suggests that this nematic-flat phase might indeed be 



unstable to thermal fluctuations. However, we have also 
shown that for D < 3, new in-plane elastic nonlinearities 
which is not included in our analysis near D = 4, be- 
come qualitatively important, and can very well lead to 
the sufficient enhancement of the bending rigidity suffi- 
cient to stabilize a nematic-flat phase of a 2D membrane. 
At the moment, this question remains unresolved. 

[22] R. Podgornik, Phys. Rev. E 52, 5170(1995). 

[23] However, topological defects in these two systems are 
qualitatively different. 

[24] As the word "phantom" suggests, we will completely ig- 
nore the effects of self-avoiding interaction. It should be 
made clear that the reason we ignore self- avoidance in 
this work is not its unimportance. On the contrary, past 
research experience in this field shows that except for 
the flat phase, self-avoiding interaction is qualitatively 
important to the long length scale properties of poly- 
merized membranes. The primary object of this paper is 
to study the global phase diagram of a liquid-crystalline 
tethered membrane. We believe that self- avoidance does 
not destroy, but rather modifies, various phases that we 
will identify in this paper. Therefore we will leave the 
technically involved study of self-avoiding interaction to 
a separate publication. 

[25] L. Radzihovsky and J. Toner, Phys. Rev. B 60 206(1999); 
Phys. Rev. Lett. 78 4414(1997). 

[26] L. Radzihovsky, A. M. Ettouhami, K. Saunders, and 
J. Toner, Phys. Rev. Lett. 87, 027001 (2001). 

[27] L. Radzihovsky and D. R. Nelson, Phys. Rev. A 44, 
3525(1991), Europhys. Lett. 16, 79(1991), Phys. Rev. 
A , 46 7474(1992); D. C. Morse, T. C. Lubensky and 
G. S. Grest, Phys. Rev. A 45, R2151(1992); D. C. Morse, 
T. C. Lubensky, Phys. Rev. A 46, 1751(1992). 

[28] L. S. Hirst and C. R. Safinya, unpublished. 

[29] Such an intrinsic nematic order parameter Q may, for 
example, characterize a spontaneous lipid tilt-order or 
alignment of liquid crystal polymers making up a nematic 
elastomer membrane. It should be intuitively clear that 
the component of such orientational order normal to the 
membrane, should have no anisotropic effects on the elas- 
ticity of the tethered membrane. Consequently Q the ne- 
matic order parameter is a well-defined reference-space 
rank-two tensor that can couple to membrane's metric 
tensor g. This contrasts qualitatively with distinct physi- 
cal realization of, for example, a tethered membrane fluc- 
tuating inside a nematic liquid crystal. In this latter case 
the appropriate nematic order parameter is a tensor in 
the embedding space that can only couple to tangent vec- 
tors of a membrane. 

[30] P G. de Gennes and J. Prost, The Physics of Liquid Crys- 
tals, Clarendon Press, Oxford, 1993. 

[31] This crumpling transition is second order at a mean field 
level and is supported by numerical simulation. RG anal- 
ysis, however, shows that it may be discontinuous. 

[32] This is mathematically identical to what's called a "min- 
imal" coupling of a charged-field and the corresponding 



12 



gauge field in gauge theories that ensures gauge invari- 
ance, with the identification of a with the gauge field, u 
the phase of the charged field and l/a^ the charge. 
[33] T. Salditt, I. Koltover, J. O. Radler, and C. R. Safinya, 
Phys. Rev. E 58, 889(1998). 



[34] This is because in two dimensions the metric tensor g is 
a 2 x 2 matrix and has only two rotationally-invariant 
independent degrees of freedom, i.e. two eigenvalues. 



